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The truncation of the perturbative series at one loop order for the mass renormalization constants remain a 
significant systematic uncertainty in the determination of heavy quark masses in lattice QCD. We present here 
a high beta Monte Carlo calculation of the two loop mass renormalization constant for clover-improved fermions 
near the charm mass in the Fermilab heavy quark formalism. A preliminary value for the charm quark mass in 
the MS scheme at two loop order is reported. 



1. INTRODUCTION 

Precision calculation of the standard model pa- 
rameters is one of the goals of lattice QCD. The 
determination of the charm quark mass, however, 
remained a difficulty as amg ~ 1. 

El-Khadra eta^j^ have shown that the Wil- 
son action has a stable heavy quark limit. At fi- 
nite lattice spacing and these heavy quark masses, 
there are different mass parameters that can be 
defined with varying lattice artifacts. By appro- 
priately tuning the input parameters and remov- 
ing the lattice artifacts by continuum limit ex- 
trapolation, one can recover the continuum mass 
definition. The agreement of the different lat- 
tice masses has been demonstrated for the charm 
quark mass by Kronfeld|^ where the matching 
to the conventional MS scheme was performed 
at the one loop level. In this work, we cal- 
culate the lattice two loop renormalization con- 
stant using large /3 Monte-Carlo techniques which 
were demonstrated for Wilson quarks in an earlier 
work|^. These constants are then used to deter- 
mine the charm quark mass in the MS scheme at 
two loop order. 

2. CHARM QUARK POLE MASS AND 
THE MS MASS 

The matching of the lattice to continuum mass 
is performed using the pole mass as an intermedi- 
ate step. In determining the pole mass from the 
lattice at finite lattice spacing, one has several 
choices for the masses that differ in their lattice 
artifacts. The rest mass (Mi) and the kinetic 



mass (M2) were used in Ref. to determine the 
MS mass at one loop order. The rest mass is de- 
rived from the spin averaged binding energy of 
the charmonium IS states and the other from the 
quark's kinetic mass. It was shown in that pa- 
per that the two quantities defined as such do in- 
deed agree in the continuum limit. In this work, 
we choose to work with the former quantity. Mi, 
mainly for its simplicity. 

The two loop relation between the pole mass 
and the MS mass has been calculated by Gray et. 
al^. The relationship in the quenched approxi- 
mation is given by, 
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where M is the pole mass. The corresponding 
expression for the pole mass determined from the 
spin averaged binding energy is M , 



M = i{M:r-a-i(^ies)*'^^ + 2a-iMl°]} 



Note that the definition of the perturbative coeffi- 
cients (M-[*') is slightly different than those given 

in Ref. (5). 

The charm quark mass in the MS scheme is 
determined at two loop order from these two ex- 
pressions. 



3. PERTURBATIVE COEFFICIENTS IN 
THE LATTICE THEORY 

Mertens et al.^ have tabulated the mass de- 
pendent one loop coefficients, m[^\ for a wide 
range of bare quark masses. We determine the 
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two loop coefficent, M[ \ using large /3 Monte 
Carlo methods H]. The procedure is essentially 
the same as described in Ref. except that 
we use a different gauge and periodic twisted 
boundary conditions. We first fix to the axial 
gauge instead of starting from the Coulomb gauge 
for these clover- improved quarks. The periodic 
twisted b.c. gives much smaller finite volume de- 
pendence than the anti-periodic b.c. used for the 
Wilson quarks. 

3.1. Simulation 

We use three values of the hopping parame- 
ter which were tuned to the spin-averaged cc ki- 
netic mass. The three values correspond to tun- 
ing performed at /3 = 5.7, 5.9 and 6.1. Mean link 
tadpole-improved tree level value for the clover 
coefficient was used in the tuning. Since we do 
not know the nonperturbative value of the crit- 
ical mass at large /3, the tadpole-improved bare 
quark mass (toq) is held fixed in the perturbative 
expansion. This gives rise to a difference in the 
treatment of the power divergences between this 
work and Ref. Q where this is treated nonpertur- 
batively. However, the knowledge of the two loop 
constant for the critical mass will allow to switch 
between the two schemes. 

The couplings and the lattice volumes used are 
similar to the ones used for the Wilson quark 
study. This includes 7 values of /3 > 9 and 7 vol- 
umes, from 4^ X 16 to 16*. The statistics is com- 
parable to the anti-periodic lattices. For /3 — 60, 
we have an extra lattice, 24*, with lower statistics 
to provide a check for the infinite volume extrap- 
olation. 

3.2. Fitting Procedure 

The quark propagators are fit to a hyperbolic 
cosine from timeslice 3 to extract the rest mass. 
Ml. These were all good fits with x^'s less than 
the degrees of freedom. The fits from timeslice 4 
and 5 gave results which were consistent but with 



K 


Mr 


m[^^ from |l 


m[^^ large /3 


0.119 


0.6318 


0.4502 


0.43(2) 


0.1227 


0.5102 


0.3195 


0.30(2) 


0.126 


0.4061 


0.1818 


0.20(2) 



Table 1 

Fit results without any constraints. 

larger errors. 

At each weak coupling, we extrapolate the rest 
masses to the infinite volume limit where we as- 
sume the following form for the finite volume cor- 
rections, 

where vi and V2 are varied independently. 

We then fit the rest masses from each weak cou- 
pling to a power series (fourth order) in av{q*)- 
The prescription for the determination of q* given 
in Ref.Q results in a very large q* especially for 
K = 0.126 since the one-loop coefficient is rather 
small. However, we note that the two loop coef- 
ficients extracted at two different reasonable q*s 
scales within the errors. We therefore use the val- 
ues of the BLM-LM combined q*s from Ref. Q. 

We summarize the result of the fits without fix- 
ing any of the coefficents in Table ^ We are able 
to reproduce the one loop coefficients within the 
statistical errors. 

3.3. Constrained Curve Fitting 

The higher order terms in the unconstrained fit 
produces coefficients that are much larger than 
the one-loop term, but also has large errors such 
that they are consistent with zero. To determine 
the two loop coefficients, we use a constrained 
fitting method using Gaussian priors |^ to reduce 
the errors. The one-loop term is constrained to 
roughly 1% of the known value. The higher or- 
der terms have been constrained with a central 
prior of with a width of 100. The fits are stable 
against the choice of the widths of the priors, but 
we note that the (regular) x^/d.o.f. for k = 0.126 
is somewhat large (> 3). The results are in Ta- 
ble 0. 



K 








0.119 

0.1227 

0.126 


0.6318 
0.5102 
0.4061 


0.452(7) 
0.320(6) 
0.185(3) 


1.3(6) 
1.4(5) 
2.1(4) 



Table 2 

Two-loop results from the constrained fits. Only 
statistical errors are quoted. 
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Figure 1. The infinite volume extrapolation of 
Ml - Mj°l at /3 = 60. Note that the point L = 24 
was not used in the fit. 



4. RESULTS/DISCUSSION 

We use the two loop results obtained above 
from large (3 Monte-Carlo to calculate the charm 
quark mass in the MS scheme at two loop or- 
der with given lattice cutoffs. The masses are 
extrapolated to the continuum limit using a? er- 
rors. The results for the one-loop, both in this 
work and Ref . , and the two loop MS mass are 
summarized in Fig. |^. Note that the difference 
in treatment of rric is significant at the one loop 
level. Our preliminary result for the two-loop 
charm quark mass in the MS scheme is 1.27(5) 
where only the statistical errors are quoted. The 
fitted two loop values are fairly stable against 
the choice of q* , truncation of the series and the 
widths of the priors. The exception is the light- 
est bare quark mass where the one loop coefficient 
is so far not reproducable without a sizeable . 
The infinite volume extrapolation is still under 
study as well as the conversion to the nonpertur- 
bative treatment of the critical mass. 

ACKNOWLEDGEMENTS 

This work was done in collaboration with 
A. Kronfeld, P. Lepage, P. Mackenzie, J. Simone 
and H. Trottier. Fermilab is operated by Univer- 
sity Research Association, Inc. for the U. S. De- 
partment of Energy. 



> 
o 




0.2 0.4 0.6 

afp_is(GeV-2) 

Figure 2. The 1-loop and 2-loop MS masses at 
finite lattice spacings. The continuum limit is 
taken with a? errors. 
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